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Abstract 

The cross sections for exclusive D^D^ and D^D^ meson pair production in peripheral nucleus - 
nucleus collisions are calculated and several differential distributions are presented. The calculation 
of the elementary 77 — t- DD cross section is done within the heavy-quark approximation and in the 
Brodsky- Lapage formalism with distribution amplitudes describing recent CLEO data on leptonic 
D^ decay. Realistic (Fourier transform of charge density) charge form factors of nuclei are used to 
generate photon flux factors. Absorption effects are discussed and quantified. The cross sections 
of a few nb are predicted for RHIC and of a few hundreds of nb for LHC with details depending 
on the approximation made in calculating elementary 77 — )■ DD cross sections. 
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I. INTRODUCTION 



The main aim of the heavy ion program at uhrarelativistic colhsion energies is concen- 
trated on the discovery and analysis of the quark-gluon plasma. High charges of the colliding 
ions give a possibility to study also peripheral processes when only a few particles are pro- 
duced. In particular, large fluxes of photons associated with the huge charges of nuclei open 
an interesting possibility to study photon-photon collisions [l[ which are difficult to study 
in e^e~ and proton-proton collisons. _ 

In our earlier works we have studied a production of [2] , p°p° [1] and recently of cc 

and hh 01 . We have demonstrated how important is the inclusion of the realistic charge form 
factors responsible for generating photon flux factors. In our analysis charge form factors 
of nuclei are calculated as Fourier transform of the realistic charge densities as measured in 
electron scattering off nuclei. 

Many years ago Brodsky and Lepage have suggested how to calculate large-angle produc- 
tion of light meson pairs (vrvr, KK) in photon- photon collisions. In the meantime several 
further studies have been preformed and some improvements have been suggested. Parallel 
those processes were searched for in several experimental studies. 

Heavy quark meson pair production was studied theoretically only in Ref. (sj where for- 
mulas have been derived in the heavy quark approximation with Dirac delta-like distribution 
amplitudes. On the other hand both lattice QCD joj and the CLEO collaboration f\ ex- 
tracted the D-meson distribution amplitude which turned out to differ considerably from 
the delta-like distribution amplitude assumed in heavy quark approximation. In the present 
studies we will use also the more realistic distribution amplitudes. 

Both LEP2 and Belle studies were not able to extract corresponding cross sections. The 
nuclear processes look potentially interesting in this respect due to the large charges of 
the colliding ions. In the present analysis we wish to calculate cross sections for exclusive 
production of D^D~ meson pairs in A A — )■ AADD reactions at RHIC (Au+Au, Wnn = 200 
GeV) and LHC (Pb-|-Pb, Wnn = 5.5 TeV). Whether the experimental analysis is possible 
requires experimental feasibility studies. 

In Fig{T]we show the basic QED mechanism of the exclusive production of DD pairs in 
the peripheral heavy-ion collisions. We consider both D^D^ and D^D^ production. 




FIG. 1: The diagram for the exclusive DD meson production. 

In the next section we concentrate on the calculation of the cross section for the 77 — DD 
reaction which, as discussed in section [Tll is used in the equivalent photon approximation in 
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the impact parameter space (b-space EPA). The nuclear cross sections are presented in the 
result section. Conclusions and summary close our paper. 

II. HEAVY CHARMED MESON PAIR PRODUCTION IN PHOTON- PHOTON 
COLLISIONS 

In this section we discuss how to calculate the elementary cross section for the 77 — ?• DD 
reaction. For the rather heavy mesons a pQCD approach seems the best method to calculate 
the cross section. In the leading order of a^, 20 Feynman diagrams are involved as shown 
in Fig|2J The diagrams can be classified into three groups. Six diagrams of first part in 
FigJ2] represent the heavier- quark pair production by two- photon collisions followed by one 
virtual gluon emission to allow the produced heavier quarks to hadronize into heavy mesons. 
The second part of the figures shows diagrams obtained by exchanging heavier quark lines 
with lighter antiquark lines. The last part consists of eight diagrams where one photon 
produces a pair of heavier quarks and the other photon produces a pair of lighter quarks. 




FIG. 2: Heavy meson pair production in photon-photon collision. The wiggly (curly) line is for a 
photon (gluon), while the solid (dotted) one is for the heavier (lighter) (anti)quark. 

In Ref.jH the distribution amplitudes of the form: 

$m(x,Q2) = /(x,Q2)l±Z^5,$^^(a;,Q2) ^ J*(^^g2)l_X^^ (2.1) 

have been assumed with the heavy-quark approximation 

/(x, Q2) = fix, Q') = 6ix - A/M) , (2.2) 

where A = M — mg. 
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The resulting pseudoscalar-pseudoscalar (PP) production amplitude Mpp(A,A') can be 
written as lal 



M'Pp{X,X') = 2- 



l-z 



2\2 



[l + z' 



elFpp{\, A') + elF^,pi\ A') - 2eQe,Gpp(A, A' 



where 



(1 - z^) 



1 — X 



X 



elHpp{\\') + ^elH'pp{\\' 



5X^(1 — xY 



M 



M 



(2.3) 



(2.4) 



Fpp(A,A') 



:i-x)[2-x{s + 2)] + -ax,-x' 



(2.5) 



X- [2- (l-x)(s + 2)]-fTA,-A' 



(2.6) 



Gpp{\ A') =[-[l + z' + {l- z^)ax,-y] - 2x{l - x(2 + s)) 
(/32-z2)-[s(l-^2)-2(3-^2)aA,_A' 



(2.7) 



Hpp{\ A') = [2 - x{s + 2)] (1 - x)/3' - z^) + ^Va,-a' + x{s + 2)aA, 



(2.8) 



H'pp{\ A') = [2 - (1 - x)(S + 2)] [x(/32 - z^) + ^Va,-a 

+ (1 -X) +X(S + 2)(TA,-A'• 



(2.9) 



In the formulas above: s 



M2 



, z = (3cos9, and (3 



1 - 



The 9 is the scattering angle between the photon and a heavy meson, the color factor Cp = 

In the e^e~ — ?■ e"'"e~7r+7r^ reaction usually cuts on 2; = cos6, where 6 is the angle in the 
photon-photon center of mass frame, are imposed. This type of cuts is rather difficult to 
realize in the AA — )■ AADD reactions where kinematics is usually incomplete. However, very 
often transverse momenta of charged pions can be easily measured with modern detectors. 
Therefore in the following we impose cuts on transverse momenta of both charmed mesons. 
In order not to loose too much of the signal we require only pt > 1 GeV. 

In Figl3] we show the energy dependence of the elementary 77 — >■ D^D^{D^D^) cross 
sections calculated according to the formalism used in Ref. j5| with the cut on meson 
transverse momenta. The cross section for D^D'^ is larger than that for D^D~ close to the 
threshold but falls faster with increasing W^^. 
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FIG. 3: The energy dependence of the 77 —5- D^D^ (black hne) and 77 — D^D^ (dashed, red 
onhne) in the heavy-quark approximation [5|]. In this calculations we have imposed an extra cut 
on the D-meson transverse momenta pt > 1 GeV. 



Above we have used the heavy quark approximation, i.e. we have assumed that light 
quark/ ant iquark carries a fixed fraction of the final meson momentum, x = ^ = ^^^^ and 
used hard process matrix elements with finite heavy quark masses. At present, it is known 
that the distribution amplitude is not of the Dirac delta type (see e.g. fio!] and references 
therein). 



Therefore in the following we use also the classical Brodsky-Lepage formalism with 
distribution amplitudes taken from the recent studies |10] . In this approach the normalized 
to unity distribution amplitude is given as 



dV 



in 



D 



x(l 



1 - Erfi 



boy 



x(l — x) 



exp 



{xm*2 + (1 - x)mf - y^) 



xil 



X] 



(2.10) 

where y = xm2 + (1 — x)m^, Erf{x) = Jq exp{—t'^)dt and ml and are quark (con- 



stituent) masses. We take the parameters in Eq. fl2.10p from [lOj: Pd — 0.8, Ad 
116 GeV, bn = 0.592 GeV-\ fo = 0.223 GeV. 

In the Brodsky-Lepage formalism the amplitude for the 77 — ?■ DD process can be written 
schematically as: 



M^^^DD = J n'^)Hx^X:,{x,y, z)^*{y, f/) dxdy . 



(2.11) 



Above Hx^\2{x,y, z) is the photon helicity dependent hard matrix element for 77 — )■ qQqQ 
and $'s are Brodsky-Lepage distribution amplitudes. The distribution amplitude $ here is 
related to the distribution amplitude in Eq. fl2.10p as: 



D 



0d(x) 



(2.12) 



In FigJH we show elementary cross section as a function of the 77 subsystem energy for 
the Brodsky-Lepage formalism. As in the heavy- quark approximation an extra cut on the 



5 



D meson transverse momentum, pt > 1 GeV, has been imposed. While the cross section for 
D~^D~ is similar in both approaches, the cross sections for D'^D^ differ considerably. 
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FIG. 4: The energy dependence of the 77 — )• D~^D~ (black line) and 77 — )• D^D^ (dashed, 
red online) in the Brodsky-Lepage formalism jTlj with distribution amplitude from 10]. In this 
calculations we have imposed an extra cut on the D-meson transverse momenta {pt > 1 GeV). 



III. NUCLEAR COLLISIONS, EQUIVALENT PHOTON APPROXIMATION 

To calculate the cross section of the nuclear process it is convenient to introduce a new 
kinematic variable: 2; = J^, where uj is the energy of the photon and the energy of the 
nucleus Ea = 'jArriproton = iMa, where Ma is the mass of the nucleus and 7 is the Lorentz 
factor. 

The total cross section can be calculated by the convolution [2|: 

a (^AA — )■ AADD; saaJ ~ J ^ ^ DD; W^^ = -y/xixiiZl) driyy (xi, X2, b) . (3.1) 

The effective photon fluxes can be expressed through the electric fields generated by the 
nuclei: 

1 1 

dra^^ (xi,X2,b) = -d^bi|E (xi, bi) |^-d^b2|E (x2, b2) |^ 

TT TT 

X 5L(b)5(^)(b-b, + b2)^^. (3.2) 

Xi X2 

The presence of the absorption factor S'^^^ (b) assures that we consider only peripheral 
collisions, when the nuclei do not undergo nuclear breakup. In the first approximation this 
can be expressed as: 

(b) = 0{h- 2Ra) = 9 (|bi - b2| - 2Ra) . (3.3) 
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Thus in the present case, we concentrate on processes with final nuclei in the ground state. 
The electric field strength can be expressed through the charge form factor of the nucleus: 

E (x, b) = / ^f''' \2^l2Mf - + ^'^i) ■ ^3.4) 

Next we can benefit from the following formal substitution: 

^ J d%\E (x, b) |2 = J d^hN (w, h)=n (u) (3.5) 

by introducting effective photon fiuxes which depend on energy of the quasireal photon co 
and the distance from the nucleus in the plane perpendicular to the nucleus motion b . 
Then, the luminosity function can be expressed in term of the photon flux factors attributed 
to each of the nuclei 



dn^^ {oou 002, b) = [e{\hi-h2\- 2Ra) N (wi, bi) N {002, ba) d^bid^bs — — . (3.6) 

J COi UJ2 



The total cross section for the AA — AADD process can be factorized into the equivalent 
photons spectra ( n {u) ) and the 77 — t- DD subprocess cross section as (see e.g.Q): 



a (^AA AADD; saa) = ^ ^ (77 ^ DD; W^^) 9 (|bi - bs] - 



2R, 



X N{co,, bi) N {U2, bs) d^bid^bs — — , (3.7) 

UJi U2 



where Wy-y = y/AuiU2 is energy in the 77 subsystem. 

a (AA AADD) = f a (-f-f ^ DD; ^/AiJ^) n (wi) n (wa) . (3.8) 

^ ' J ^ ' UJy UJ2 

Additionally, we deflne Y = \ {iid + Vd)) rapidity of the outgoing meson system which is 
produced in the photon-photon collision. Performing the following transformations: 

= = ^e-^ , (3.9) 

^"^^"^ '-diy,,dy, (3.10) 



Ul U2 W-y-y 

d{ui,U2) 



dujiduj2 — )• dW^^dF where 
formula (13. 7p can be rewritten as: 



d{W,„Y) 



(3.11) 



cr 



(^AA AADD; Saa) = j ^ (ll ^ DD; W^^) 9 (|bi - bsl - 2R^ 



xN{uj^M)Niuj2,h2)T^d%id^h2dWyydY . (3.12) 

Finally, the nuclear cross section can be expressed as the flve-fold integral: 

a [aA ^ AADD; saa) = J a {-f-f ^ DD; W^^) 9 (|bi - ba] - 2Ra) 

xN{coiM)N{u2,h2)2nb^dbmd\dby^dWyydY , (3.13) 
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TABLE I: Total cross section for the exclusive DD production for RHIC (g/y/y = 200 GeV) and 
LHC (sjvAf ~ 5500 GeV) calculated with distribution amplitudes from 0] and 



Process 


CTtot 




BCS94 DA 


ZH2007 DA 


AuAu AuAuD^D^ 


12.44 nb 


0.08 fib 


AuAu AuAuD^D^ 


21.2 nb 


5.06 nb 


PbPb PhPhD+D- 


1.68 fib 


1.92 fib 


PbPb PbPbD^D^ 


4.09 fib 


0.28 fib 



where bx = {bi^ + &2a)/2, by = {biy + &2j/)/2 and bm = bi — 62 have been introduced. This 
formula is used to calculate the total cross section for the AA — )■ AADD reaction as well as 
the distributions in the impact parameter b = b^, the meson invariant mass W^^ = Mjjf, 
and the meson pair rapidity Y{DD). 



IV. RESULTS FOR THE NUCLEAR COLLISIONS 

In this section we present results for gold-gold coUisions at RHIC {Wn^ = 200 GeV) and 
lead-lead collisions at LHC (Wnn = 5.5 TeV). In all the calculations presented here we use 
realistic charge form factor being a Fourier transform of the realistic charge densities. Both 
D^D~ and D^D^ channels are discussed. Below we present quantities which can be easily 
calculated in the b-space EPA discussed in the previous section. In Table I we have collected 
numerical values of the cross sections for considered exclusive processes. The cross sections 
are small but the reactions seem measurable especially at LHC. 

Before we go to the presentation of more interesting quantities we wish to show a control 
plot. In Figjniwe show distributions in the auxiliary variables b^ and by defined in Section UTTl 
These distributions are peaked at bx,by = and strongly decrease with increasing \bx\ or \by\. 
This figure demonstrates that the chosen range of integration over the auxiliary variables 
bx and by is sufficient. Please note that the range of integration for the LHC energy is 5 
times larger than that for the RHIC energy. To speed up calculation we have performed 
integration in only two quadrants in the [b^, by) space and used relevant symmetry relations. 

In the heavy-ion photon-photon processes particles can be produced when nuclei fly far 
away apart. This is particularly spectacular for light leptons (see e.g. j2|). In FigE] we 
show such a distribution for the pseudoscalar DD production. The maximum of the cross 
section is for the impact parameter b when colliding nuclei almost touch each other and the 
cross section decrease for larger b. The decrease is much sharper for RHIC than for LHC. 
Compared to light leptons j2| here it is easier to achieve a convergence in b. 

Now we come to the distributions which can be directly measured in experiments. In 
FigJ7]we show the distribution in the D^D~ and D^D^ pair rapidity being Y ^ + 
The visible irregularities at larger |y| are caused by the oscillating nuclear form factor. 
The larger meson pair rapidity the larger four-momentum squared of the exchanged photon 
(which is the argument of the charge form factor). The irregularities correspond to the four- 
momentum squared when the charge form factor changes its sign, i.e. when |Fp is close to 
zero. A more elaborate discussion on the role of the nuclear form factor can be found in 
Ref.g. 
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FIG. 5: The distributions in the auxihary quantities bx and by defined in section Hill for D^D 
and Z)°Z)° for the BCS94 DA [5] (black) and for the ZH2007 DA ^lOj] (red) for RHIC and LHC 
energies. 



In FiglH] we show distribution in invariant mass of D and D. We predict steep fall-off 
of the disribution as a function of the invariant mass, even steeper than for the elementary 
cross section shown in FigJ31 This figure shows that in experiments only a region of small 
invariant masses close to the threshold could be investigated. 



V. CONCLUSIONS 



We have calculated for the first time in the literature total and differential cross sections 
for exclusive production of pseudoscalar DD meson pairs in the AA — )■ AADD reaction 
assuming that the reaction is driven by the 77 — )■ DD subproceess. The elementary cross 
sections were calculated in the heavy - quark approach as well as in the Brodsky- Lapage 
formalism with distribution amplitude describing recent CLEO data on leptonic D^ decay. 
Rather small cross sections have been found. The cross section for exclusive DD produc- 
tion is much smaller than the cross section for the exclusive or semi-exclusive production 
of cc calculated recently. In our calculations absorption effects were included in the impact 
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FIG. 6: Distribution in impact parameter for D^D and D^D^ for the BCS94 DA [5] (black) and 
for the ZH2007 DA [10,] (red) for RHIC and LHC energies. 
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FIG. 7: The cross section as a function of the meson pair rapidity for D^D and D^D^ for the 
BCS94 DA 01 (black) and for the ZH2007 DA ^ (red) for RHIC and LHC energies. 



parameter Equivalent Photon Approximation. The meson pairs are produce preferentially 
when the nuclei almost touch each other. The cross section strongly depends on the ap- 
proximation made in the calculation. The dominant contribution to the cross section comes 
from the region of very small DD invariant masses. 

Whether the process can be measured requires further dedicated studies including de- 
tector simulations. We belive that our evaluation will be a useful starting point for such 
future studies. Since the cross section strongly depends on theoretical details experimental 
verification would be very helpful. 
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